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ABSTRACT 

Aims. By means of optical Integral Field Spectroscopy observations, we aim to disentangle and characterize the starburst component 
in the Blue Compact Dwarf Galaxy Mrk 1418. In particular we propose to study the stellar and ionized gas morphology, to investigate 
the ionization mechanism(s) acting in the interstellar medium, to derive the physical parameters and abundances of the ionized gas. 
Methods. Integral Field Spectroscopy observations of Mrk 1418 were carried out with the Potsdam Multi-Aperture Spectrophotometer 
(PMAS) at the 3.5 m telescope at Calar Alto Observatory. The central 16" x 16" (1.14 x 1.14 kpc^ at the distance of Mrk 1418) were 
mapped with a spatial sampling of 1"; we took data in the 3590-6996 A spectral range, with a linear dispersion of 3.2 A per pixel. 
The seeing was about r.'5. From these data we built maps of the most prominent emission lines, namely [On], H/?, [Oiii], Ho-, [Nii] 
and [Sn] as well as of several continuum bands, plus maps of the main line ratios: [O iii]/IIy6, [Niil/Hcf, [Sii]/Hq', and KafH/j, and 
derived the physical parameters and gaseous metal abundances of the different star-forming regions detected in the field of view. 
Results. Mrk 1418 shows a distorted morphology both in the continuum and in the ionized gas maps; the current star- formation 
episode is taking place in five knots, distributed around the nucleus of the galaxy. The interstellar medium surrounding these knots 
is photo-ionized by stars, with no clear evidence for other excitation mechanisms. The galaxy displays an inhomogeneous dust 
distribution, with the high Ha/H^S ratio in the central areas indicating a large amount of dust. The oxygen abundances derived for 
the individual star-forming knots are very similar, suggesting that the ionized interstellar medium is chemically homogeneous in O/H 
over spatial scales of hundreds of parsecs. This abundance (Z a 0.4Zo from the empirical calibrations) places Mrk 1418 among the 
high metallicity BCD group. 

Conclusions. These findings show the advantages of IFS when investigating the properties of such complex objects as BCDs, with 
an asymmetric star forming component. Only a bidimensional mapping of their central regions allows to approach such questions as 
the star formation processes in BCDs, the star-forming history of the individual starburst knots, or the abundance gradients. 

Key words. galaxies:dwarf ~ galaxies:starburst ~ galaxies: stellar populations -galaxies:individual (Mrkl418) -ISM: HII regions 



K*" 1. Introduction 



Blue Compact Dwarf (BCD) galaxies have a unique potential to 
address many central issues in contemporary galaxy research. 
Chemically unevolved nearby star-forming (SF) dwarfs, such 
as BCDs, are in fact an important link to the early Universe 
and the epoch of galaxy formation, as they are regarded as 
the local counterparts of the distant subgalactic units (build- 
ing bl ocks) from which larger systems are created a t high red- 
shifts dKauffmann et alJll993l:lLowenthal et al.lll997l) : the study 
of these systems may provide important insights into the star- 
formation process of distant galaxies. Moreover, even though 
most BCDs are not genuinely young galaxies, their metal defi- 
ciency make them useful objects to set constraints on the primor- 
dial "^He abundance, and to monitor the synthesis a nd dispersal 
of heavy ele ments in a near l y pristine env ironment {Pagel et aL 
ll992; Maseg osa et aLlll994l:Eotov et al.i ri997; Kunth & Ostlin 
[2000 ). BCDs are also ideal laboratories for the study of the star- 



burst phenomenon: being smaller and less massive than normal 
galaxies, they cannot sustain a spiral density wave and do not 
suffer from disk instabilities, which considerably simplifies the 
study of the star-formation process. Besides, the radiation emit- 
ted by their SF regions is less affected by the stellar continuum 
than in giant spiral galaxies, allowing for more accurate studies 
of element abundances. 

In spite of the great effort spent on the study of BCDs dur- 
ing the last two decades, many fundamental questions such as 
the mechanisms responsible for the ignition of the starburst, 
their evolutionary status, or their star formation histories, are 
still far from being well understood. Detailed spectrophotomet- 
ric analyses of individual nearby BCDs, providing the informa- 
tion needed to disentangle young and old stellar populations, as 
well as the gas and the dust components, are crucial to further 
invest igate these topics (Papader os et al.lll998l: IGil de Paz et al 
2000t (Cairos etalj|2002; .Guseva et all 12003 



120071) . 



bll3; Cairos et al 



* Based on observations obtained at the German-Spanish 
Astronomical Center, Calar Alto, operated by the Max-PIanck- 
Institut fur Astronomic Heidelberg jointly with the Spanish National 
Commission for Astronomy. 



Yet this ki nd of analysis is scarcely found in the literature 
(see lCairos et al. 2007, and references therein). This fact is prob- 
ably due to the large amount of time needed to obtain the re- 
quired data using traditional observing techniques, i.e. the com- 
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bination of broad and narrow-band imaging and long-slit spec- 
troscopy; this can easily lead to observing times of several nights 
per object. Moreover, these observations typically suffer from 
varying atmospherics conditions, which can render the combi- 
nation of the various data sets a complicated and tedious under- 
taking. 

The relatively new observational method of Integral Field 
Spectroscopy (IFS) allows to overcome these problems. Each 
single exposure with an Integral Field Unit (IFU) provides at 
the same time both spatial and spectral information, making IFU 
observations an order of magnitude more efficient than tradi- 
tional observing techniques. Furthermore, IFS provides simul- 
taneous spectra for all spatial resolution elements, under the 
same instrumental and atmospheric conditions, resulting in a 
homogeneou s dataset. As several recent studie s have shown 
dlzotoy et al.1 120061: iGarcfa-Lorenzo et al.l l2008t iKehrig et all 
120081 IVanzi etalJl2008h Integral Field Spectrographs are ideal 
instruments to study such small and yet complex objects as 
BCDs. 

We have undertaken a comprehensive spectrophotometric 
analysis of a large and representative sample of BCDs by means 
of IFS. This sample, which includes w 40 objects, covers the 
whole range in luminosity (from Mb = -21 to Mg - -14) 
and morphologies (nE, iE, IC, and il, following the classifica- 
tion scheme introduced bv lLoose & Thuaiil (Il986l) ) that we find 
among BCD galaxies. 

The dataset resulting from this study will help answer some 
of the most crucial open questions in BCD research: to effec- 
tively disentangle the young and old stellar populations, to set 
constraints on the age of the galaxies, to study the kinematics 
of the gas and stars, to look for signatures of merger or interac- 
tion episodes, and, finally, to determine their star-forming his- 
tories. Besides, these data will provide the essential template to 
understand and interpret the results for intermediate and high-z 
star-forming galaxies. 

This is the first of a series of papers in which the outcomes 
from this project will be presented. In this paper we carry out a 
pilot study of the galaxy Mrk 1418, with the intention of showing 
the great potential of IFS when applied to the BCDs field. Mrk 
1418 is part of a sample of ten objects which were observed with 
the Potsdam Multi-Aperture Spectrophotometer (PMAS) at the 
3.5 m telescope at Calar Alto; results for the remaining galaxies 
observed with PMAS will be presented in a forthcoming paper 

Mrk 1418 is a low luminosity (Mb - -16. 96) BCD which 
accord ing to the morphological classification bv lLoose & ThuaiJ 
(Il986l) . belongs to the most common BCD class, the iE 
BC Ds. It has bee n included in t he BCDs sampl e s ana lyzed 
by iDoubher etaP ([l997l IT999I) . iGil de Paz et all (l2003h and 
iGil de Paz & Madord ( l2005h as w ell as in the recently pu blished 
"Atlas of Markarian galaxies" bv lPetrosian et alj (|2007[) . Its ba- 
sic data are shown in Table [1] Its proximity (distance of 14.6 
Mpc) and the asymmetric, clumpy starburst morphology make 
this object an excellent test-bench to illustrate the need of bidi- 
mensional spectroscopy to analyze BCD galaxies. Its current 
SF activity takes place in a number of knots placed on a regu- 
lar a nd elliptical en velope (B - R ^ 1.6; Doublier et al. 199^ 
iGil de Paz & Madord l2005l) . Ivan Zee et alj (120011) found, from 
high spatial resolution H i synthesis observations, that the neutral 
gas in Mrk 1418 extends to approximately two times the optical 
diameter of the galaxy (D25 - 48 arcsec, Df|j/D25 = 2.3, see 
Table 1 in van Zee et al. 2001) and displays kinematical pecu- 
liarities reminiscent of tidal tails. Althoug h the optical profi le is 
better fitted by a de Vaucouleurs law (.Doubfier et alJ[T997l) . the 



HI radio observations presented in lvan Zee et alj ( 1200 ll) indicate 
that the galaxy is supported by rotation. 

The paper is structured as follows: in Sect. |2] we describe 
the observations, the data reduction process and the method em- 
ployed to derive the two-dimensional maps. In Sect.[3]we present 
the main results of this work, that is, the emission line and con- 
tinuum flux, emission line ratio, and velocity maps, as well as 
the results derived from the analysis of the integrated spectra 
of the selected galaxy regions. These results are discussed and 
summarized in Sect. |4] and Sect.|5] respectively. 



2. The data 

2.1. Observations 

Mrk 1418 was observed in 2007 March with the Potsdam Multi- 
Aperture Spectrophotometer (PMAS) at the 3.5 m telescope 
at the Observatorio Astronomico Hispano Aleman, Calar Alto 
(CAHA). PMAS is an integral field spectrograph, with a lens ar- 
ray of 16 X 16 square elements, connected to a bundle of fiber 
optics, whose 256 fibers are re-arranged to form a pseudoslit 
in the focal plane of the spectrograph (see Roth et al. 2005 and 
iKelz etal.ll2006i for more details about the instrument). In the 
configuration we used, each lenslet covers 1 " x 1 " in the sky, 
thus providing a total field of view of 16" x 16" (1.14 x 1.14 
kpc^ at Mrk 1418 distance, with a spatial sampling of 71 pc/") . 

A grating with 300 grooves per mm was used during the 
observations, in combination with a 2048 x 4096 pixel SITe 
ST002A CCD detector. This set-up provides a spectral range of 
3590-6996 A, with a linear dispersion of 3.2 A per pixel (the 
CCD was binned 2 x 2 in the spatial and dispersion directions). 
The actual spectral resolution, obtained by measuring the width 
of the arc lines in comparison spectra, was about 6.8 A FWHM, 
corresponding to ~ 300 km s"' at Ha and ~ 400 km s ' at O iii. 

We observed a total of 4800 s on the galaxy (four expo- 
sures of 1200 s each); additional sky frames were obtained by 
off'setting the telescope several arcmin away from the target. 
Calibration frames were taken before and after the galaxy ex- 
posures. These calibrations consist of spectra of emission line 
lamps (HgNe lamp), which are required to perform the wave- 
length calibration, and spectra of a continuum lamp, necessary 
to locate the individual spectra on the CCD and to perform the 
flat-fielding correction. As usual, bias and sky-flats exposures 
were taken at the beginning and at the end of the night. The 
spectrophotometric standard stars BD+75325 and BD+332642 
were observed for flux calibration. The seeing, as indicated by 
the CAHA seeing monitor, was about r.'5 FWHM. 



2.2. Data reduction 

The data have been processed using standard IRAlQ tasks. The 
reduction procedure includes bias subtraction and image trim- 
ming, tracing and extraction, wavelength and distortion calibra- 
tion, flat-fielding, combination of the individual frames, sky- 
subtraction and flux calibration. 

The first step in the data reduction was the bias subtraction. 
All the bias exposures were averaged to obtain a master bias, 
which was then subtracted from the rest of the frames. 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Next, apertures (that is, the spectra produced by each fiber) 
were defined on well exposed continuum frames using the task 
apall; the task first finds the centroid of the aperture at a fixed 
position on the dispersion axis, and then asks for the width of 
the extraction window, which we set at 6.4 pixels (the best com- 
promise between maximizing the flux from each spectrum and 
avoiding contamination by nearby fibers). The apertures were 
then traced by fitting a polynomial to the centroid along the dis- 
persion axis (because of the field distortion introduced by the op- 
tics of the system, the apertures show a noticeable curvature, de- 
pending on their spatial position). A fifth degree Legendre poly- 
nomial was found to provide good fits, with a typical RMS of 
about 0.01 pixels. 

Once the apertures were defined and traced in the continuum 
frames, we used apall again to extract them in all the images. 
The extraction consists on summing the pixels along the spatial 
direction into a final one-dimensional spectrum. The final prod- 
uct is the so-called "collapse" spectrum: an image MxN, where 
M is the number of pixels in the dispersion direction, and is 
the number of spaxels (256 for PMAS). 

Afterward we performed the wavelength calibration and the 
dispersion correction. In order to calibrate in wavelength we 
used the tasks identify and reidentify: i) first, in the comparison 
spectra (arc) we identified several emission features of known 
wavelength in a reference fiber; ii) second, a fifth degree poly- 
nomial was fitted across the dispersion direction, resulting in a 
RMS of about one hundredth of A; iii) next, with reidentify we 
identified the emission lines in the remaining fibers of the arc 
frame, using the selected one as reference. The polynomial fit 
derived in this way was used to calibrate the rest of the images 
(galaxy, standard stars and sky frames) using the task dispcor. 

The data were corrected for throughput (that is, spaxel-to- 
spaxel overall sensitivity variations) and response (detector pixel 
non-uniformity) using the task msrespId. The individual galaxy 
spectra were then corrected for atmospheric extincti on, adopting 
the "su mmer extinction coefficients" published bv .Sanchez et al.l 
(l2007h . and combined using the task evicombine. 

The next step was the sky subtraction. Sky frames (5 minutes 
exposures) were processed in the same way as galaxy spectra. A 
one-dimensional sky spectrum was produced by averaging, with 
a sigma-clipping algorithm, the signal along the spatial direction. 
The flux of the three to four brightest sky lines was measured in 
both the final sky spectrum and in the final galaxy spectrum to 
determine the scale factor by which to multiply the sky spectrum 
before subtracting it from the galaxy spectrum. Because the rela- 
tive intensity of different sky lines vary noticeably on short time 
scales ([Patat 2003), it is very difficult to find a scaling factor that 
applies equally well to all the sky lines and the sky continuum, 
and some fine tuning is required. The final scale factor was found 
by trial and error, aiming at minimizing overall residuals in the 
sky lines (especially those close to galaxy emission or absorp- 
tion lines) even if it leaves large residuals in bright sky lines not 
affecting any interesting spectral features. 

Finally, the flux calibration was performed. The spectra of 
the two spectrophotometric standards observed the same night as 
Mrk 1418 were reduced in the same way as the galaxy frame, ex- 
cept that the sky spectrum to be subtracted out was computed by 
median-combining the outermost fibers. The one-dimensional 
spectra of the standard stars were obtained by summing all the 
fibers within a radius of about 2 FWHM (typically 3 to 4 arcsec- 
onds) from the fiber with highest signal. The tasks standard and 
sensfunc were used to derive the sensitivity curves, after com- 
bining together the data for the different spectrophotometric stars 
observed in the same night. 



By comparing the sensitivity curves obtained for all the spec- 
trophotometric s tandards observed throughout the three nights 
of the observing run, we estimate that the relative uncertainty on 
the calibration factor is generally equal or less than 2%, except 
blueward of 4000 A, where the curves show a marked change of 
slope and the uncertainty increases up to about 8%. 

2.3. Emission Line Fits 

In order to measure the relevant parameters of the emission lines 
(center, flux and width), they were fitted by a single Gaussian. 
The fit was carried out by using the minimization algorithm 
implemented by C. B. Markwardt in the mpfitexpr IDL librarjB 

The continuum (typically 30-50 A on both sides) was fitted 
by a straight line. Lines in a doublet were fitted forcing them to 
have the same redshift and width. 

Criteria such as flux, error on flux, velocity, and width were 
used to do an automatic assessment of whether to accept or re- 
ject a fit. For instance, lines with too small (< 2.5 A) or too 
large (> 5.0 A) widths were flagged as rejected, as well as lines 
with relative flux errors ^ 10% — the exact limits depending on 
the specific line. Such criteria were complemented by a visual 
inspection of all fits, which led to override, in a few cases (typi- 
cally 5-10 spaxels), the automated criteria decision. For instance, 
sometimes the program fitted large spikes for very faint of even 
completely absent emission lines; in other cases, a particularly 
noisy continuum or bad pixels decreased the computed S /N be- 
low the acceptance threshold, while the actual fit was clearly 
good enough. 

For the H/3 line, where absorption wings were present, we 
simultaneously fitted two gaussians, representing the absorption 
and the emission components. To improve the reliability and sta- 
bility of the fit, an additional constraint was applied to the equiv- 
alent width of the absorption component, by allowing it to vary 
only between and 5 A. 

2.4. Maps construction 

The emission line fit procedure gives for each line parameter (for 
instance flux) a table with the fiber ID number, the measured 
value and the acceptance/rejection flag. This table was then used 
to produce a 2D map, by using an IRAF script that takes ad- 
vantage of the fact that PMAS fibers are arranged in a regular 
16 X 16 matrix. The script automatically converts ADU counts 
into flux (erg s"' cm"^ A"'), by multiplying by the wavelength 
dependent conversion factor taken from the sensitivity curve of 
the flux calibration. 

Two continuum maps were obtained by summing the flux 
within specific wavelength intervals (4500 to 4700 Aand 6000 to 
6200 A), selected so as to avoid emission lines or strong resid- 
uals from the sky spectrum subtraction. Line ratios maps were 
simply derived by dividing the corresponding flux maps. 

3. Results 

3. 1 . Morplioiogy of tlie ionized gas and stars 

To study the galaxy morphology we use flux maps of the bright- 
est emission lines and the two continuum bands above. The 
emission line maps trace the interstellar gas, ionized by massive 



^ |http://cow.physics.wisc.edu/~craigni/idI/idl.htmll 
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Fig. 1. Left panel: B-band image of Mrk 1418 with isocontours, spaced 0.5 mag apart, overlaid. Right panel: continuum subtracted 
Ha frame, with the 16 x 1 6 arcsec^ PMAS field o f view overlaid and the main star-forming knots labeled. North is up, east to the 
left. Both images ai-e from lGil de Paz et all (l2003h . and pubHcly available in the NED; the PSF is 2'.'6 FWHM. 



Table 1. Basic data for Mrk 1418 



Parameter 



Value 



RA (2000) 
DEC (2000) 
Vhciio (kms"') 
D (Mpc) 
Ab (mag)" 
triB (mag)*" 
inR (mag)'' 
Mb (mag)" 
FQla) (ergs 
mj (mag)"" 
m/f (mag)"" 
iuk, (magf 
Mm (Mef 



cm 



09:40:27 
+48:20:16 
773 ± 11 
14.6+ 1.0 
0.098 

13.86 + 0.03 
12.61 +0.19 
-16.96 
95 + 16 X 10- 
12.08 + 0.03 
11.46 + 0.04 
11.41+0.05 
2.30 X 10** 



Notes: RA, DEC, heliocentric velocity, Vheiio, and distance, D (Mpc), 
all taken from NED (http://nedwww.ipac.caltech.edu/ 1. Distance cal- 
culated using a Hubble constant of 73 km s"' Mpc"', and taking into 
account the influence of the Virgo Cluster, the Great Attractor and the 
Shapley supercluster. 



Absorption coefficient in the B band, from lSchlegel et al] ( |1998|) . 
B- and i?-band integrated magnitudes, their Icr un certainty and in- 
tegrat ed Hot flux, corrected for Galactic extinction jGil de Paz et alj 
|2003l) . 

Absolute magnitude in the B band, computed from the above tabu- 
lated B integrated magnitude and distances. 

Total magnitudes from 2M ASS. 

Neutral hydrogen mass Mm from Ivan Zee et all ( |2001|) . which 
adopted a distance of 10.8 Mpc. 



stars, while the continuum maps show the stellar light distribu- 
tion. For comparison. Figure [T] displa ys the published B broad- 
band and Ha images of Mrk 1418 (Gil de Paz etalJl2003) . with 
the field of view covered by the PMAS observations overplotted. 
There is good qualitative agreement between the reconstructed 
maps (see Figure© and the direct images. 



Figure|2]displays the emission line flux maps in [O ii] /13727, 
HyS, [Oiii] /t5007. Ha, [Nii] A65M, and [Sii] ^^6717, 6731. 
The galaxy shows the same morphology in all emission lines: 
five SF knots are distributed in a roughly circular pattern around 
the optical center of the galaxy. The emission peak is located 
in knot 1 in all the lines (see Figure [TJ, in the western part of 
the galaxy. This large knot seems to be connected with knot 4, 
placed » 10" to the northeast, with a bridge-like structure. 

In the continuum (see bottom panels in Figure |2]), Mrk 1418 
shows an inner kidney-shaped feature, with a peak located ap- 
proximately at the center of the outer isophotes. The position 
of this central maximum in the continuum frames does not de- 
pend on what specific "emission line free" spectral range is used, 
hence we can safely define this continuum peak as the optical 
nucleus of the galaxy. Smaller condensations are detected north 
and southeast of the galaxy nucleus. 

The structure of the ionized gas, as traced by the emission 
lines, looks quite different from the stellar light distribution. At 
the position of the continuum peak, there is a dearth of gas emis- 
sion, indicating that the star formation in the central regions has 
already ceased. The position of the knots detected in emission 
does not coincide with any continuum peak (knot 2 is close, but 
not coincident with the continuum blob seen in the southeast cor- 
ner). Knots 3 and 4 seem to be co-located with a shallow, local 
blue maximum in the "blue-red" continuum map, while knots 1 
and 2 are slightly displaced from the blue continuum peaks. 

3.2. Ionization: mechanism and structure 

To obtain information about the ionization structure of the 
ionized nebulae and the nature of the ionization sources 
acting in Mrk 1418, we built emission line ratio maps in 
[Oiii] A5QQl/n/3, [Nii] A65MfHa and [Sii] AA67n, 6731/Hq'. 
Although [O i] /16300 is detectable in many spaxels, it unfortu- 
nately falls very close to the sky line [O i] /16300, which prevents 
us from obtaining precise measurement of its flux and, hence, of 
die [Oi] /16300/Ha ratio. 

Using the information derived from the above mentioned line 
ratios it is possible to distinguish among the different ionization 
mechanisms acting in a gaseous nebula, namely, photoioniza- 
tion by radiation from OB stars, shock-wave heating and pho- 
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Fig. 2. Emission line maps for [On] .t3727, H/3, [Om] A5QQ1, Ha, [Nii] A65U and [S ii] AA6in, 6731, continuum maps in the 
intervals 4500^700 A (blue) and 6000-6200 A (red), and the continuum color map. Axis units are arcseconds; north is up, east 
to the left. Maps are in logarithmic scale (in units of 10"'^ erg s"' cm"^), except the "blue-red" continuum color map which is in 
magnitudes (with arbitrary zero point). The black dot marks the nucleus of the galaxy, defined as the peak in the continuum frames. 
The Ha map also shows the outline of the SF knots and of the two "interknot" regions, A and B (see text for details). The isophotes 
of the "red" continuum map are spaced 0.2 dex apart. 



toionization by a power law con tinuum source dBaldwin et al.l 
1 198 It IVeilleux & Osterbrockll 198 7). The [Oiii] ^15007/^8 emis- 
sion line ratio is an excitation parameter indicator and provides 
information about the available fraction of hard ionizing pho- 
tons of the ionizing star cluster embedded in the nebula; there- 
fore a large [O iii] A5QQ7/H/3 ratio indicates a highly ionized re- 
gion. [Nil] /16584 and [S ii] AA61\1 , 6731 are low ionization 
lines, usually weak in Hii regions: large [Nii] /16584/Ha' and 
[Sii] AA61\1 , 6731/Hq' ratios usually indicate the presence of 
ionizing mechanisms dif ferent from photoioniz ation, e.g. shock 
waves and/or AGNs (e.g. lShuU & McKedl 19791) . 



Figure [3] displays the maps 
[Nil] A65Ufi\a and [Sii] ^^6717, 



of [Oiii] ^5007/H/3, 
6731/Hq'. In order to 



allow a quick comparison with the position of the ionization 
sources. Via contours have been overlaid. These three maps 
show basically the same structure. As expected in regions 
photoionized by stars, they trace the star-formation activity, 
and the peaks in [Om] /15007/Hj6 coincide with the minima 
in [Nil] /16584/Hq' and [Sii] ^^6717, 6731/Ha. Only knot 4 
displays a different behavior: it is a moderate minimum in 
[Nil] ^6584/Ha and [Sii] ^^6717, 6731/Hq' and also shows 
low values in [Om] ^5007/Hy3. The [Om] ASQQllWp map also 
displays a peak at about (-2", -2"), the same position as the 
peak in the extinction map, in an area that does not coincide 
with any emission knot. 
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Fig.3. Line ratios maps: [Oiii] ASOQl/Hfi, [Nii] A65M/Ha, [Sn] AA6717, 6731/Hq', Ha/Hfi and [Sii] /l6717/[Sii] ^6731. The 
last panel shows the Ha equivalent width map. Ha contours (spaced 0.2 dex apart) have been overlaid in all the maps. Axis units 
are arcseconds; north is up, east to the left. All maps are in logarithmic scale, except Ha/H/S. The black dot marks the nucleus of the 
galaxy. 



3.3. Extinction 

The extinction map was built from the Ha/H/? ratio. No reliable 
fits could be obtained, for the individual spaxels, for H6 and Hy, 
as they are much weaker and are also strongly affected by under- 
lying stellar absorption. 

The Ha/HyS ratio map is shown in Figure[3] Mrk 1418 shows 
an inhomogeneous extinction pattern, which peaks in the cen- 
tral, inter-knot regions, displaced from the continuum maximum 
(x! 3" to the southeast). There, Ha/H/S varies between 4 and 5.5, 
suggesting strong extinction by dust with a reddening E{B - V) 
of up to 0.62 mag, while in the surrounding SF knots the HafH/3 
ratio is < 3, that is E{B - V) < 0.04 mag. Indeed, an extinction 
patch or lane located between the emissi on line knots has been 
alrea dy observed in several BCDs (e.g. Garcia-Lorenzo et alj 
I2OO&: Vanzi et al. 20 08; Kehrig et al. 200&) , and could be related 
to the sweeping out of the ISM by stellar winds and SNs. In the 
continuum "blue-red" color map, the position of the extinction 
peak is redder than the knots, but no dust lane or patch is clearly 
seen; changes in stellar populations as well as the specific geom- 
etry of the dust distribution may account for this seeming incon- 
gruity. 

This result stresses the importance of performing a bidimen- 
sional correction of the interstellar extinction: just assuming a 
single, spatially constant value for the extinction can lead to 
large errors in the derivation of fluxes and magnitudes in the 
different regions of the galaxy. As an example, if we used the 
extinction value for knot 1, E(B - V) ~ 0.07, to correct the 
observed colors in the nuclear region (where the extinction is 



E(B -V)~ 0.40) we would obtain a color B - V about 0.3 mag 
redder than its actual value. 



3.4. Electron density distribution 

We also produced the map of the ratio [S n] /l6717/[Sii] /16731, 
an electron number density sensitive diagnostic for the range 
100-10.000 cm"^ (FigureEll. We find that the density is almost 
constant, with values close to the low density limit in almost the 
whole field of view. Only in a few fibers in the central, inter-knot 
region (the same one where the peaks of the [O iii] ASQQl/Hfi ra- 
tio and of the extinction are located) the density is higher. 

3.5. Integrated spectroscopy of selected galaxy regions 

As the next step in our spectroscopic analysis, we used the de- 
rived flux and continuum maps to isolate the regions of interest 
in the central parts of the galaxy. They include the five bright- 
est SF knots detected in the narrow-band maps and the nucleus, 
defined as the continuum peak (see Figure|2l). 

We identified the fibers corresponding to each region, and 
summed their spectra together to create the final spectrum. The 
regions were defined as those spaxels with an Ha flux larger than 
3.2 (in logarithm), a threshold chosen as the best compromise be- 
tween enclosing a sufficient number of spaxels for each knot, and 
avoiding that the knots overlap or abut each other A few spaxels 
above this threshold but whose inclusion would have stretched 
the shape of a knot were excluded, while knot 5 was augmented 
by two additional spaxels below the threshold. 
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3.5.1. Reddening corrected line fluxes 
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Fig. 4. Spectra of the 5 knots, of the nucleus, and the total spec- 
trum obtained by summing all the fibers. Spectra are shown in 
logarithmic scale and are offset for clarity. Small tickmarks are 
spaced 0.2 dex apart. 



We also defined two interknot regions: A, where the Ha flux 
has a local minimum, and B, where the [S ii] AA&lll, 6731/Ha 
ratio is highest (see Figure|2]). 

In this way we produced higher signal-to-noise spectra (in 
comparison to the spectra of the individual fibers), better suited 
to derive physical parameters and chemical abundances. The 
area of these regions ranges between ~ 0.02 and 0.10 kpcl We 
also computed the total spectrum of Mrk 1418 by summing all 
the 256 spaxel spectra (that is, over an area of 1.3 kpc^). 

Fig. |4] displays the spectra for the six selected regions and 
the whole mapped field of view. The spectra of the five SF re- 
gions show similar characteristics: all of them are dominated by 
strong emission features, superimposed on a blue continuum. 
Strong [On] Xilll , 1^8, [Oiii] A5Q01 , Ha, [Nn] /16584 and 
[Sii] AA6in, 6731, are visible in emission. Absorption wings 
are detected in H(5, Hy and HjS, for knots 3, 4, and 5. Although 
[Oi] /16300 is detected in the spectra of all the five knots, no 
reliable flux measurements could be obtained, as the profile of 
this line is contaminated by the residuals from the subtraction 
of the [Oi] /16300 sky-line. The spectrum of the nucleus shows 
strong absorption features, indicating an important contribution 
of older stars; H6, Hy and higher order Balmer lines (HI 1, HIO, 
H9, H8, H7) are visible only in absorption. 



Fluxes and equivalent widths of the emission lines were mea- 
sured by fitting the line profile with a Gaussian. Reddening- 
corrected intensity ratios and equivalent widths for the different 
spatial regions are listed in Table |2l The correction from under- 
lying stellar absorption was done in two different ways, as ex- 
plained below. 

In knots 4 and 5, and in the nuclear and the integrated spec- 
tra, the Hy6 Balmer line shows clear absorption wings, and we 
fitted simultaneously an absorption and an emission component. 

In the remaining cases, we proceeded as follows. We first 
adopted an initial estimate for the absorption equivalent width, 
EWabs, corrected the measured fluxes, and computed the extinc- 
tion coefficient C(Hy6) through a least-square fit to the Balmer 
decrement given by the equation: 



F(A) 1(A) 



Fim I(Hf3) 



(1) 



where 



F{H/3) 



ized to H^; 



is the line flux corrected for absorption and normal- 
-M^ is the theoretical value for case B recombina- 



tion, from iBrocklehursJ (Il97l l). and /(^ is the reddening curve 
normalized to H/3 which we took from lWhitfordI (Il958l) . 

We then varied the value of EWabs, until we find the one that 
provides the best match (e.g. the minimum scatter in the above 
relation) between the corrected and the theoretical line ratios. 
This is done separately for each knot. 

We note that this correction accounts for the fraction of 
emission line flux that is masked by the underlying absorption, 
and not necessarily represents the actual absorption equivalent 
width. To illustrate this point, let us consider two (extreme) ex- 
amples: (1) the width of the emission line and that of the absorp- 
tion line are very similar; in this case the absorption component 
is completely hidden underneath the stronger emission compo- 
nent, and the measured emission flux will be the algebraic sum 
of the fluxes of the two components (the flux of the absorption 
line being negative); (2) the absorption line is very broad and 
shallow; in this case, the emission line flux measured by setting 
the continuum level on both sides of the "base" of the emission 
line profile will only slightly underestimate the actual flux, and 
the correction will be much smaller that the actual equivalent 
width of the absorption component. 

We compared our measurements for the integrated and 
the nuclear spectra with the spatially integrated (within a 
30" X 50" area) and nuclear ( 2" 5 x 2"5) values published by 
iMoustakas & KennicuttI (l2006l) . The results are presented in 
Fig. |5] keeping into account the different sizes of the apertures 
used, and a likely small spatial offset in the location of the nu- 
clear aperture, the comparison looks fair. 

3.5.2. Line ratios, physical parameters and chemical 
abundances 

Table [3] lists the most relevant line ratios, physical parameters 
and chemical abundances for the selected regions in Mrk 1418. 
In order to calculate the errors associated with the physical prop- 
erties and chemical abundances, we propagated the emission- 
line intensity errors listed in Table |2] 

The physical conditions (electron density and temperature) 
and ionic element abundances were derived from the reddening- 
corrected emission line fluxes followi ng the 5 -level atom five I 
program in the IRAF n ebular package jPe Robertis et al.ll 19871 ; 
IShaw & Dufou3[T995h . 
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Fig. 6. Optical emission-line diagnostic diagram for the differ- 
ent SF knots and the two interknot regions (A and B). The 
curves separate Seyfert galaxies, LINERS and Hn regions- 
like objects; solid lines are the empirical boundaries from 
IVeiUeux & Osterbrockl (Il987h. while dotted lines are the theo- 
retical boundaries from lKewley et alJ (l200lh . 



Fig. 5. Comparison between the integrated and nuclear fluxes 
and equiv alent widths from this paper (TP ), and the values pub- 
Ushed by iMoustakas & Kennicutl (l2006h (M&K), which used 
apertures of 30" x 50" and 2'.'5 x 2"5 respectively. For this com- 
parison, our measurements have only been corrected for galactic 
extinction and, for the Balmer lines, for underlying absorption 
(see text for details). Filled circles represent the integrated spec- 
tra, open circles correspond to the nuclear spectra. The line at 
^ 4350 A is Hy. 



Electron densities were calculated from the ratio of 
[Sii] A6in/A613l emission lines. The electron tempera- 
ture of [Om] (re[Oiii]), representative of the high-excitation 
zone of the ionized gas, was derived from the [Om] 
A4363/(/l4959-H^5007) ratio for knots 1 and 2, and for the 
summed spectrum. In order to calculate re[Oii] we used the 
relation between [O n] and [Oiii] electron temperatures from 
iPilvugin et al.l (l2006l) . We assumed the approximation re[Sii] 
a! re[Nn] ^ T^[Oii] for the calculation of S+ and N+ abun- 
dances, since no auroral line could be measured in the low exci- 
tation zone. We adopted re[Om] for the calculation of O^^. As 
the nebular He ii /14686 A line was not detected in any of the knot 
spectra, we can assume that the contribution from highly ionized 
species like O^^ is negligible. Therefore, the total oxygen and ni- 
trogen abundances were obtained as: O/H = (O^ H-0^^)/H^ and 
N/O = N+/0+. 

To obtain oxygen abundances in those knots in which 
[Om] could not be measured (knots 3, 4, 5, and the 
nucleus), we apphed the c ommo nly used strong-line 
method from iPettini & Pagell ( |2004|) . These authors re- 
vised the N2 {= log[Nii] A65^4/Ua) and C>3A^2 (= 
log{([Om]/l5007/H/3)/([Nii]^6584/Ha'))) indices, using 137 



extragalactic Hii regions. The uncertainties in the derived 
metallicities from these two indices are ~ 0.38 and 0.25 dex 
respectively. 

For knots 1 and 2, the estimated oxygen abundances from 
empirical parameters are consistent, within the uncertainties, 
with those derived using the electron temperature. 

The position of the five SF regions in the diagnos- 
tic diagram: [Om] .15007/H^ vs [Nii] ^6584/Ha' and 
[Sii] AA67n, 6731/Hq', is shown in Fig. |6] The em- 
pirical boundaries between the different zones (from 
I Veilleux & Osterbrocki i 19871). as well as the theoretical 
boundaries proposed by iKewlev et al.l (1200 ll) . are also plotted. 
In both diagrams, all the SF knots fall within the area occupied 
by Hii regions. 

3.6. Kinematics of tlie ionized gas 

We obtained the velocity field by fitting a single Gaussian to the 
strong Ha and [Om] A5Q01 emission lines. 

Results are displayed in Fig. |7] No regular pattern is visible 
in these maps, whereas some rotation is expected in the inner- 
most 10", where the velocity c hanges from 760 to 800 km/s in 
the radio data by lvan Zee et alj jToOl) (see their Fig. 12). Given 
the spectral resolution of our data, ~ 300 km s ' at Ha and ~ 400 
km s"' at O iii, velocity structures can not be reliably measured. 
An analysis of the kinematics of Mrk 1418 in the optical has to 
await data taken with higher spectral resolution. 

4. Discussion 

Results based on PMAS IFS observations of the BCD Mrk 1418 
are presented. The mapped area covers the central starburst re- 
gion of the galaxy (16 x 16 arcsec^, equal to 1.14 x 1.14 kpc^; 
see Fig.lTJ. 
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Fig. 7. Velocity field of the ionized gas in the central region of the galaxy for the [O iii] /15007 and the Ha lines. Axis units are 
arcseconds; north is up, east to the left. 



The emission line and continuum maps reveal quite an in- 
triguing structure. Both gas and stars display an iiTegular mor- 
phology; interestingly, none of the peaks in the ionized gas 
maps coincides with a maximum in the continuum. Small spa- 
tial offsets between continuum and emission line peaks are 
common in compact starburst ga laxies, as reported for inst ance 
by [H unter & Thronson ( 1995), Maiz-Apellaniz et al.1 (Il998l) and 
iLagos et al.l (|2007), and are likely related to the release of kinetic 
energy by massive stars and supernova explosions. 

In the emission lines, at least five major knots are resolved; 
the two biggest ones, knot 1 and knot 2, have approximate diam- 
eters of 400 pc and 300 pc respectively. These sizes are com- 
parable to those typical of giant H ii regions in M 33 (e.g., 
NGC 604'Sandage & Tammani] 11974 and in the LMC (e.g., 
30Dor;IWalborn 1991). 

In the continuum maps a bright nucleus is located in the cen- 
tral part of the field of view, in a place with reduced Ha emission 
(see Figure |2|l. In the continuum the galaxy shows a distorted, 
kidney-shaped morphology, with boxy isocontours, which can 
be seen with better clarity in the B-band image in Fig. [T] and 
in the red continuum image in Fig. |2] (see contours). Such boxy 
morpholog y, though usually interpreted as a sign of interactions 
or mergers jPoublier et al.|[l997lll999l) . can just be the result of 
the superposition of the starburst on an older, more regular stellar 
component. 

The current SF activity is taking place in the five major line- 
emitting knots (Fig. 12]), while the lack of significant gas emission 
in the nuclear region suggests the presence of a dominating in- 
termediate age populatio n. The underlying hos t galaxy detected 
in br oad-band images dPoublier et alj 119971; i Gil de Paz et alj 



2003h shows colors (B - R ~ 1.5 - 1.7; |Doublier et alj[l997t 



Gil de Paz & Madorg|2005l) indicative of a population several 



Gyrs old. 

The ages of the selected regions have been estimated by 
comparing their Ha and H/3 equivalent with the predictions 
from Starburst 99 (Leithereret al. 1999) evolutionary synthe- 
sis models. We chose those models with Z - 0.4 Zo, the closest 
value to the metallicity derived from the emission-line fluxes — 
a reasonable approximation to the metallicity of a young pop- 
ulation — and a Salpeter initial mass function (IMF) between 
1 and 100 Mq. We found that we can reproduce the measured 
Ha and Hy6 equivalent widths for the five SF knots with an in- 
stantaneous burst of star formation and ages ranging between 5 



and 7 Myr: knots 1 and 2 seem to be the youngest, with equiv- 
alent widths consistent with ages between 5 and 6 Myr, while 
for knots 3, 4 and 5, we found values slightly larger, from 6 to 
7 Myr. The equivalent widths in the nucleus are consistent with 
ages of about 10 Myr. 

We should keep in mind that these values are indeed upper 
limits to the ages of the knots. In fact, as the continuum flux 
increases due to the contribution from the older stars, the mea- 
sured equivalent widths decrease and this in turn results into 
larger derived a ges for th e knots. As shown in previous works 
dCairos et alJl2(j02L l2007l) . this effect is not negligible in BCDs, 
though difficult to quantify: it rests on the model adopted to de- 
scribe the underlying component and can strongly vary with the 
position of the knots across the galaxy. For instance, the correc- 
tions found for the different knots in the galaxy Mrk 35 range 
from 10% to 90% dCairos et alj|2007l) . By assuming a contribu- 
tion of the host galaxy of 50% to the continuum in the brightest 
knots in Mrk 1418, the computed age will be about 1 Myr less. 

The Interstellar Medium (ISM) around the SF knots is 
excited by the ionizing photons coming from the central 
stars/clusters, with no clear evidence for other excitation mech- 
anisms. In the "interknot" area B a moderate [Oiii] ASOQl/Hfi 
ratio together with a large [Sn] AA6in, 6731 /Ha ratio (see 
Fig. |6l) suggest that some mechanisms other than photoioniza- 
tion contribute to the gas excitation (most probably shocks; e.g. 
Shull & McKee 1979; Schwartz et al. 1993). 

On the other hand, nebular He ii /14686 A emission has been 
detected (in one just spaxel) in a location southwest of the galac- 
tic center (near knot 2). This line offers a powerful diagnostic 
of the ionizing fluxes and stellar population in starburst regions 
but, unfortunately, here its S /N is too poor for a reliable mea- 
surements of its FWHM and flux. 

The interstellar extinction has an overall irregular distribu- 
tion, with a ridge that crosses the galaxy in a north-south di- 
rection between knots 2 and 5, consistent with material swept 
out by stellar winds and SNe. Its position and shape are in 
broad agreement with a dust-lane visible in the SDSS frames. 
This finding is in line with recent publications ([H unt et al. 200 Vl 
all I200I IVanzi & Sauvagd llool ICdros et all 120071) 



ICairos et 
that I 



^i & Sauvagd 12004 ICdros et all 1200"/ 
tnat questioned the prev ious belief that BCDs have little or no 
dust. iLerov et al.l ( l2005l) report Mrk 1418 as a nondetection in 
their CO survey. 
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In agreement with s everal other works on BCDs 
(e.g. 'Izotov et al.' '1997; 'Vi'lchez & Iglesias^PMain^ Il998t 
[Papaderos et al. 2006; Kehrig et al. 2008), we did not find any 
significant chemical abundance variation among the diff'erent 
knots in the galaxy, which points to a common chemical 
evolution scenario in all of them. Also, the nitrogen-to-oxygen 
abundance ratio values found are consistent with those obtained 
for BCDs and Hii galaxies CVflchez & Iglesias-Paramo 2003; 
iPerez-Montero & DfazllIOOSh . 

We found for the total spectrum a lower value of direct 
oxygen abundance and a higher electron temperature than for 
knots 1 and 2 (see Table |3]l. While this finding may reflect a 
genuine difference between the oxygen abundance in the total 
spectrum (kpc-sized aperture) and in the individual H n regions 
(with sizes of the order of 100 pc), it may also be due to the rela- 
tively large measurement uncertainties for the weak [O in] A4363 
emission line. 

The origin of the starburst episodes in SF dwarf galaxies is 
far from being well understood. One of the current scenarios 
is the stocha stic self-propagating star formation (SSPSF) model 
proposed bv^ Gerola et al.l (11980'), in which the starburst episode 
in dwarfs is explained as statistical fluctuations in the SFR. In 
Mrk 1418 both the ages of the starburst and the morphology 
speak against this scenario; the different SF knots have ages too 
similar to be causally related: with the typical distances among 
knots (400-800 pc), the star formation will need times of the 
order of 40-80 Myr to propagate through shock traveling at a 
typical speed of ~ 10 km s 

An alternative scenario, which is gaining popularity, is inter- 
actions and mergers. Recent searches have shown that a substan- 
tial fraction of SF dwarfs possess low-mass c ompanion galaxies, 
most of them gas-rich b ut faint in the optical (Tavlor et al.l ll993i 
I1995L 119961; iNoeskee tal. 2001; Pustilnik et al. 2001), while 
studies focused on individual objects have shown that interac- 
tions indeed play a substantial role in the evolution of these sys- 
tems (Johnso n et alj|2004l;l0sthn et al.ll2004tlBravo-Alfaro et al 



2003). 



20041 i2006t ICumming et aTl |2008| - but see IVaduvescu et al 



These observational results have been recently backed up 
from the theoretical point of view. The numerical simulations 
presented in Bekki (2008) can satisfactory explain the physical 
properties of BCDs (e.g. blue compact core, redder low surface 
brightness component, presence of compact young stellar clus- 
ters, low metallicity and extended HI gas disk) as the result of 
mergings between dwarfs with larger fraction of gas and ex- 
tended gas disks. 

In Mrk 1418, the asymmetric spatial distribution of the SF 
knots, their different sizes, the bridge-like structure connecting 
knot 1 and knot 4, and the filamentary structures visible in the 
NED Ho- image, as well as the boxy, kidney-shape morphology 
of the continuum, are all suggestive of a system that has recently 
experi enced some interac tion events. Kinematic results from H i 
maps dvan Zee et al.ll200lb also agree with a merger/interaction 
scenario: the extended gaseous distribution show kinematical pe- 
culiarities, reminiscent of tidal tails. The H i velocity field dis- 
plays a complex behavior in the inner regions, but an overall 
clear rotation pattern. 

The low velocity resolution of our PMAS data hampers the 
derivations of a reliable velocity distribution. Deeper, higher 
spectral resolution IFU observations, aimed at measuring the gas 
velocity dispersion as well as the stellar kinematics, are essential 
to put tighter constraints on the galaxy dynamics. 



5. Conclusions 

The central 16" x 16" region of the BCD galaxy Mrk 1418 has 
been mapped by using the PMAS spectrograph. Integral Field 
Spectroscopy allowed us to simultaneously study the morphol- 
ogy, ionization mechanism, extinction, chemical abundances and 
ionized gas kinematics of the galaxy. From our analysis we high- 
light the following results: 

1. Mrk 1418 displays very different morphologies in the contin- 
uum and in the emission lines. While in the continuum the 
central region has a distorted, kidney-like shape and boxy 
isocontours, with the peak located approximately at the cen- 
ter of the outer isophotes, the emission line maps exhibit five 
SF regions irregularly distributed around the optical center 

2. The excitation ratios [Oiii]/Hy6 and [Nii]/Hq' are typical of 
Hii regions in the whole mapped area, while in the inter- 
knots zones (region B) the [S ii]/Ha ratios have values > 0.4, 
possibly indicating the presence of shocks. 

3. The galaxy displays an inhomogeneous dust distribution, 
with a peak value Ha/Hfi = 5.5 in the region between knot 2 
and knot 5, about 3" southeast of the continuum peak; this 
dust patch has an area of about 8n". 

4. For each SF knot we derived line fluxes, interstellar extinc- 
tion, the most relevant line ratios and metal abundances. We 
find that: i) all the knots show Hii-like ionization (e.g. star 
formation); ii) the density does not vary significantly from 
knot t o knot, and has val ues typical of giant extragalactic ob- 
jects (Mc Call et al.n i985); iii) all knots have similar metal 
abundance, with no abundance variations on scales of hun- 
dreds of parsec. 

This pilot paper shows that IFS is an essential tool to study 
BCD galaxies, as, in just one shot, it allows us to map the gas and 
stars morphology, determine the physical parameters and derive 
chemical abundances. 

Therefore, IFU studies of large BCD samples will enable 
us to gain insights into the main issues concerning this class 
of galaxies, namely, to disentangle and characterize their stel- 
lar populations, to elaborate their SF histories and to investigate 
the mechanisms that trigger their star-formation activity. 
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Tkble 2. Reddening corrected line ratios, normalized to ^S'*\ for the knots identified in Mrk 1418. 
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82.14 ± 0.93 


3.43 ± 0.30 


62.70 ± 1.03 


3.10 ± 0.31 


52.25 ± 1.81 


6.00 ± 0.56 


117.2 ± 3.53 


5.24 ± 0.57 


65.37 ± 3.65 


10.45 ± 1.38 


23.53 ± 0.78 


5.44 ± 0.43 


43.28 ± 0.25 
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7.68 ± 1.62 
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7.31 ± 0.32 


0.37 ± 0.07 


5.31 ± 0.46 


0.26 ± 0.08 
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0.46 ± 0.02 
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1.00 
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[Oin] 


1.24 ± 0.04 


38.81 ± 0.29 


1.16 ±0.03 


41.07 ± 0.44 


1.12 ± 0.05 


19.89 ± 0.41 


0.61 ± 0.04 


11.54 ± 0.48 


1.14 ±0.08 


15.78 ± 0.66 


0.83 ± 0.09 


2.61 ± 0.18 


1.06 ±0.03 


11.01 ±0.06 


5007 


[Om] 


3.62 ±0.10 


113.9 ± 0.46 


3.45 ±0.10 


122.70 ± 0.51 


3.19 ±0.12 


55.06 ± 0.49 


1.78 ± 0.09 


34.54 ± 0.54 


3.29 ±0.18 


46.81 ± 0.90 


2.34 ± 0.20 


7.61 ± 0.20 


3.16 ± 0.09 


33.48 ± 0.09 
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0.07 ± 0.01 
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2.48 ± 0.26 
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1.88 ±0.35 
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1.60 ± 0.49 
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[Nn] 
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8.15 ±0.44 
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2.65 ±0.34 
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0.10 ±0.03 


2.37 ± 0.62 


0.12 ±0.03 


0.80 ±0.21 


0.11 ±0.01 


1.95 ±0.05 
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Hff 


2.87 ±0.10 


137.6 + 0.39 


2.88 ± 0.14 


154.1 ± 1.35 


2.90 ± 0.20 


73.72 ± 0.79 


2.87 ± 0.15 


84.12± 1.41 


2.87 ±0.16 


69.03 ± 1.48 


2.86 ± 0.34 


19.63 ± 0.27 


2.86 ± 0.07 


49.87 ± 0.10 


6584 


[Niil 


0.26 ±0.01 


12.50 + 0.23 


0.29 ± 0.02 


15.22 ±0.43 


0.22 ± 0.02 


5.68 ± 0.40 


0.31 ±0.03 


9.14 ±0.61 


0.19 ±0.03 


4.56 ±0.61 


0.33 ± 0.05 


2.25 ± 0.20 


0.27 ±0.01 


4.63 ± 0.06 


6717 


[Sii| 


0.26 ±0.01 


12.50 + 0.28 


0.31+0.02 


17.05 ±0.44 


0.31 ±0.03 


8.21 ±0.44 


0.42 ± 0.03 


12.93 ±0.60 


0.29 ± 0.03 


7.30 ± 0.66 


0.58 ±0.08 


4.21 ±0.23 


0.40 ±0.01 


8.49 ± 0.05 


6731 


[Sn] 


0.21 ±0.01 


10.22 ± 0.25 


0.23 ± 0.02 


12.69 ± 0.43 


0.20 ± 0.02 


5.43 ± 0.41 


0.32 ± 0.02 


9.87 ± 0.57 


0.20 ± 0.03 


5.17 ±0.68 


0.42 ± 0.06 


3.02 ±0.21 


0.30 ±0.01 


5.52 ± 0.07 






F(HfS) = 


324.6 ± 22.7 


F(m = 


166.3 ± 15.2 


F(m = 


67.4 ± 11.3 


f(H8) = 


84.6 ± 6.0 


F{HfS)-. 


= 42.4 ± 3.9 


F(m = 


106.3 ± 30.4 


F(m = 


1871.9 ±62.6 






C(m = 


0.111 ±0.032 


cms) = 0.073 ± 0.035 


C(H8) = 0.123 ±0.064 


cms) = 0.267 ± 0.025 


C(m = 


0.361 ± 0.030 


am = 1 


B.625± 0.105 


C(m = 


0.316 ±0.011 
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= 3.0A 
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, = 2.1A 



(*) Notes. - Reddening-corrected line fluxes, normalized to FQifi) = 1. Equivalent widths of Balmer lines are corrected for underlying steUar absorption. The reddening coefHcients, CQJfi), the value of the absorption correction, EWabs, and the 
reddening-corrected H/3 flux, f (HyS)(xlO"" ergs cm"^ s~') are listed for each region. In the nuclear region, only Hfi and Hfi could be measured, so we were unable to compute EWabs and C(Hfi) simultaneously; to determine the latter we assumed 
EWabs = 3.0. The quoted uncertainties include measurements, flux-calibration and reddening coefHcient errors. 



Table 3. Physical parameters and abundances 



Parameter 


knot 1 


knot 2 


knot 3 


knot 4 


knot 5 


Nucleus 


Integrated 


, / [0 III] ,15007 \ 


0.56 ± 


0.02 


U.j4 ± U.Uz 


0.50 ± 0.03 


A O C 1 A A1 

0.z5 ± 0.03 


r\ CO 1 f\ f\A 

0.52 ± 0.04 


0.37 ± 0.05 


A CA 1 A AO 

0.50 ± 0.02 


1 /[Nil] ,16584 \ 

^ \ Ho- / 


-1.04 = 


t0.02 


-1.00 ± 0.02 


-1.12 + 0.05 


-0.97 ± 0.05 


-1.18 ±0.09 


-0.94 ± 0.08 


-1.03 ±0.02 


Jq„[ [S II] ,1,167 17. 6731 \ 


-0.79 = 


t0.02 


-0.73 + 0.02 


-0.75 + 0.03 


-0.59 ± 0.04 


-0.77 ± 0.07 


-0.46 ± 0.06 


-0.61 ±0.02 




190 


< inn 

=i lUU 


< inn 


1 ns 

iUO 


< inn 

=i iUU 


< inn 

=i lUU 
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T (\0 III) flO'* K~) 


1.18± 


0.07 


1.22 ± 0.07 










1 67 + 30 


re([Oiii])(io'*K) 


1.27 ± 


0.10 


1.33 ±0.10 










1.96 ±0.42 


12 + \og{OIH) - (r,) 


8.09 ± 


0.06 


8.03 ± 0.06 










7.71 ±0.13 


12 + log(0/H) - {N2) 


8.26 


8.28 


8.23 


8.30 


8.20 


8.31 


8.27 


12 + log(0/H) - (03^2) 


8.22 


8.24 


8.21 


8.34 


8.19 


8.31 


8.24 


12 + log(S+/H+) 


5.87 ± 


0.05 


5.90 ± 0.05 










5.80 ±0.10 


12 + log(N+/H+) 


6.52 ± 


0.05 


6.59 ± 0.05 










6.26 ±0.11 


log(N/0) 


-1.29 = 


tO.lO 


-1.18 + 0.10 










-1.27 ±0.21 



Notes: ([O ii]) derived from the relation: re([0 ll]) = 0.72 x re([0 ill]) + 0.26 found bv lPilvugin stsil (l2006h . 

Te ([Om]): electron temperature measured from [Oiii] /14363 

12 + log(0/H) - (Je): direct 0/H abundance derived from re[Oiii] 

12 + log(0/H) - (A'2): O/H derived from the N2 index (Pettini & Pagel 2004'); the associated uncertainty is ±0.38 

12 + log(0/H) - (03^2): O/H derived from the 03A'2 index tPettini & Pagelii2004) : the associated uncertainty is ±0.25. 



